A numerical experimental study of inverse preconditioning for the parallel iterative solution to 3D finite element flow equations

Integration of the subsurface flow equation by finite elements (FE) in space and finite differences (FD) in time requires the repeated solution to sparse symmetric positive definite systems of linear equations. Iterative techniques based on preconditioned conjugate gradients (PCG) are one of the most attractive tool to solve the problem on sequential computers. A present challenge is to make PCG attractive in a parallel computing environment as well. To this aim a key factor is the development of an efficient parallel preconditioner. FSAI (factorized sparse approximate inverse) and enlarged FSAI relying on the approximate inverse of the coefficient matrix appears to be a most promising parallel preconditioner. In the present paper PCG using FSAI, diagonal and pARMS (parallel algebraic recursive multilevel solvers) preconditioners is implemented on the IBM SP4/512 and CLX/768 supercomputers with up to 32 processors to solve underground flow problems of a large size. The results show that FSAI may allow for a parallel relative efficiency larger than 50% on the largest problems with p=32 processors. Moreover, FSAI turns out to be significantly less expensive and more robust than pARMS. Finally, it is shown that for p in the upper range may be much improved if PCG–FSAI is implemented on CLX.     

Estimates for Elements of Inverse Matrices for a Class of Operators with Matrices of Special Structure

In this paper, we consider questions related to the structure of inverse matrices of linear bounded operators acting in infinite-dimensional complex Banach spaces. We obtain specific estimates of elements of inverse matrices for bounded operators whose matrices have a special structure. Matrices are introduced as special operator-valued functions on an index set. The matrix structure is described by the behavior of the given function on elements of a special partition of the index set. The method used for deriving the estimates is based on an analysis of Fourier series of strongly continuous periodic functions.     

增强型延迟反馈法控制低维混沌系统的解析研究

基于时间延迟反馈控制混沌系统的方法,提出一种增强型控制方案,并利用分析延迟系统产生Hopf分支条件的方法,给出这种方案控制低维连续自治混沌系统时,在达到控制目标的条件下,控制参数的一般解析关系.将这一方案和分析方法应用到两个混沌模型中,结果表明:采用修正的方案可以明显地改善控制混沌的效果和质量;解析分析的结果与实际数值计算的结果一致. 关键词： 延迟反馈 混沌控制 Hopf分支     

~ Transform Demonstration of Dark Soliton Solutions Found by Inverse Scattering

One of the basic problems about the inverse scattering transform for solving a completely integrable nonlinear evolutions equation is to demonstrate that the Jost solutions obtained from the inverse scattering equations of Cauchy integral satisfy the Lax equations. Such a basic problem still exists in the procedure of deriving the dark soliton solutions of the NLS equation in normal dispersion with non-vanishing boundary conditions through the inverse scattering transform. In this paper, a pair of Jost solutions with same analytic properties are composed to be a 2 × 2 matrix and then another pair are introduced to be its right inverse confirmed by the Liouville theorem. As they are both 2 × 2 matrices, the right inverse should be the left inverse too, based upon which it is not difficult to show that these Jost solutions satisfy both the first and second Lax equations. As a result of compatibility condition, the dark soliton solutions definitely satisfy the NLS equation in normal dispersion with non-vanishing boundary conditions.     

Complex Behaviors of a Simple Traffic Model

In this paper, we propose a modified traffic model in which a single car moves through a sequence of traffic lights controlled by a step function instead of a sine function. In contrast to the previous work [Phys. Rev. E 70 (2004)016107], we have investigated in detail the dependence of the behavior on four parameters, ω, α, η, and a1, and given three kinds of bifurcation diagrams, which show three kinds of complex behaviors. We have found that in this model there are chaotic and complex periodic motions, as well as special singularities. We have also analyzed the characteristic of the complex period motion and the essential feature of the singularity.     

Chaotic-to-ordered state transition of cathode-sheath instabilities in DC glow discharge plasmas

Transition from chaotic to ordered state has been observed during the initial stage of a discharge in a cylindrical DC glow discharge plasma. Initially it shows a chaotic behavior but increasing the discharge voltage changes the characteristics of the discharge glow and shows a period subtraction of order 7 period → 5 period → 3 period → 1 period, i.e. the system goes to single mode through odd cycle subtraction. On further increasing the discharge voltage, the system goes through period doubling, like 1 period → 2 period → 4 period. On further increasing the voltage, the system goes to stable state through two period subtraction, like 4 period → 2 period → stable.     

LOCAL AND GLOBAL HOPF BIFURCATIONS FOR A PREDATOR-PREY SYSTEM WITH TWO DELAYS

In this paper, the Leslie predator-prey system with two delays is studied. The stability of the positive equilibrium is discussed by analyzing the associated characteristic transcendental equation. The direction and stability of the bifurcating periodic solutions are determined by applying the center manifold theorem and normal form theory. The conditions to guarantee the global existence of periodic solutions are given.     

A Stability Estimate for a Solution of the Problem of Determining the Dielectric Permittivity

We consider the problem of determining the dielectric permittivity for a nonconducting and nonmagnetic medium. As information we take the traces of the tangential components of the electromagnetic field on the lateral surface of a cylindrical domain. These traces correspond to a solution to some direct problem for the Maxwell system. The impulse source of the current flux lies outside the domain in which the coefficient is sought. The main result of the article is a stability estimate for a solution to the inverse problem in question.     

Gamow Vectors in Exactly Solvable Models

Through two examples: the Friedrichs model and a particular case of central potential scattering, we illustrate the way of constructing Gamow vectors.     

Stability and bifurcation in a delayed predator-prey system with Beddington-DeAngelis functional response

We consider a delayed predator-prey system with Beddington-DeAngelis functional response. The stability of the interior equilibrium will be studied by analyzing the associated characteristic transcendental equation. By choosing the delay τ as a bifurcation parameter, we show that Hopf bifurcation can occur as the delay τ crosses some critical values. The direction and stability of the Hopf bifurcation are investigated by following the procedure of deriving normal form given by Faria and Magalhães. An example is given and numerical simulations are performed to illustrate the obtained results.